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Abstract. Gravitational collapse of a class of spherically sym- 
metric stars are investigated. We quantise the geometries describ- 
ing the gravitational collapse by a deformation quantisation pro- 
cedure. This gives rise to noncommutative spacetimes with gravi- 
tational collapse. 
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1. Introduction 

Gravitational collapse [l]-[5] is one of the most dramatic phenom- 
ena in the universe. When the pressure is not sufficient to balance 
the gravitational attraction inside a star, the star undergoes sudden 
gravitational collapse possibly accompanied by a supernova explosion, 
reducing to a super dense object such as a neutron star or black hole. 

Theoretical investigations predicting gravitational collapse were car- 
ried out in the early 1930s in the ground breaking work of Chan- 
drasekhar [HI E]- In 1939, Oppenheimer and Snyder [8J investigated 
the collapse process of ideal spherically symmetric stars equipped with 
the Tolman metric [9] -dr^ + e^^^^'^UR^ + e^(^'^)(rf^2 ^ sin^^rf?/^^)^ 
When the energy-momentum of an ideal star is assumed to be given 
by perfect fluids, Tolman's metric allows the case of dust which has 
zero pressure. In the dust case, Oppenheimer and Snyder solved the 
Einstein field equations by further assuming that the energy density 
is constant. They showed that stars above the Tolman- Oppenheimer- 
Volkoff mass limit [SI ID] (approximately three solar masses) would col- 
lapse into black holes for reasons given by Chandrasekhar. The work 
of Oppenheimer and Snyder also marked the beginning of the modern 
theory of black holes. 

The purpose of the present paper is to investigate gravitational col- 
lapse in noncommutative general relativity. We shall work within 
the framework of the noncommutative Moyal geometry developed in 
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[TT| [T2| [13] . Gravitational collapse is in principle understood classically, 
however, quantum effects will become important at the final stage of 
the collapse, especially if a star collapses into a black hole. As there 
still lacks a consist theory of quantum gravity, it is useful to incorpo- 
rate some quantum effects into gravity by deforming general relativity. 
Much effort has been made in this area in recent years, resulting in 
several tentative proposals for noncommutative general relativity |14] - 
HZ], HB US US], [IB]-iI]. Much work has been done to investigate 
noncommutative corrections to black holes, and we refer to|12j, [23j- 
[29J for details. 

On a cautionary note, we should mention that the widely cited papers 
[I6| [T7] were shown in |22] to yield results entirely different from the low 
energy limit of string theory. If the cause of this is not mathematical 
imprecisions in [161 HI]! then it may be an indication that there are 
flaws in the rationale of these papers, as the low energy limit of string 
theory for an appropriate choice of vacuum is physically realistic. 

In the present paper, we first investigate static interior solutions of a 
class of spherically symmetric stars equipped with the Misner-Zapolsky 
metric [34j (see equation (12.11) ) and have energy-momentum given by 
perfect fluids (see equation f l2.2p ). The energy densities of the stars 
are taken to be decreasing functions instead of (physically unrealistic) 
constant functions. New solutions of the Einstein field equations are 
obtained, and their singularities signalling gravitational collapse are 
discussed. We feel that these results are interesting by themselves even 
from the point of view of classical gravity. 

We then quantise the metrics obtained by the deformation quanti- 
sation procedure developed in [HI |T21 113] • This gives rise to noncom- 
mutative spacetimes with static interior singularity. 

Finally, we quantise the dust solutions studied by Oppenheimer and 
Snyder [8]. It gives rise to the noncommutative dynamical gravita- 
tional collapse. As far as we are aware, gravitational collapse in the 
noncommutative setting has not been studied before. 

2. Spherically symmetric stars 

Suppose a spherically symmetric star is equipped with the following 
metric 

g = -e^^^'^-^e + e2^(*'^)rfr2 + r\de^ + sin^ Od^lj^), (2.1) 
and its energy-momentum is given by perfect fluids with 

T^u = {e+p)U^U,+pg^^, (2.2) 
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where e is the energy density, p the pressure and the vector U is the 4- 
velocity of the fluid elements. This metric was studied by Misner and 
Zapolsky [M] for neutron star models. Its physical properties differ 
considerably from those of the Tolman metric [9]. In particular, it 
rules out dust solutions (with p = but e 7^ 0). We shall show that, 
for general spherical symmetric metric (12. ip . if the energy density is 
independent of the time t, then the metric must be static. We shall 
also study the cases that energy densities are given by step functions 
and quadratic decay functions. 

We mention that in some textbooks, e.g. [21 El E]; the metric (12. ip 
was discussed in the static case with functions a, (3 depending on r 
only. The appearance of a curvature singularity was interpreted as an 
indication of gravitational collapse. 

We choose the co-moving coordinates and take the four-velocity U 
to be pointing in the timelike direction, that is, Ui = U2 = = 0. 
We also normalise U by setting If^Ui, = —1. Now the components of 
energy-momentum tensor in the coordinates {t, r, 6, ip) are given by 



T 



2 

r p 



sin^ I 



Denote the Einstein tensor by G^i, = R^^ — ^g^v The Einstein field 
equations give that [4] 



1 2drl3\ 1 
e I ^ ^ = -e, 



1^2 lyt J /^'2i 



^/^jLi ly J /y ^ 

Gl =e-'^3l = 0, 



r 



- (^d^f3 + {dtf3Y - dtaM') = p. 

As a consequence of the Bianchi identity, we have the Tolman- 
Oppenheimer-Volkoff (TOV) equation [H [TOl l9] 

{6 + p)dra + drP = 0. (2.3) 

We now show that the metric (12. ip does not allow any dust solution: 
the vanishing of the pressure p implies the vanishing of the energy 
density e. In fact, (12. 3p implies either e = or dra = Oifp = 0. 
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In the latter case the equation gives (3 = 0, thus the Gl equation 
gives e = 0. By the Birkhoff theorem, (12 .ip must be Schwarzschild in 
this case. This is different from the Tolman metric, for which the dust 
solution was constructed by Oppenheimer and Snyder [8]. 

It is straightforward that the Gg and G^ equations are consequence 
of those of G*, G;, G[ and the TOV equation Q (see, e.g, [21 El [5] for 
the proof of this fact in the static case). In fact, the Gl, G^ equations 
give 

' 2r ^2 2r' ' 

2r ^ 2 2r' 
From the (j[ equation, we obtain 

/3 = /3(r). 

This implies that 

e = €{r). (2.4) 
From the TOV equation (12.31) . we obtain 

-e^^(^ I ^ 3p 1 rV 



4 4 4 2r2 4 
A' 2r' 



Substituting these into the equation, we find that the left hand side 
is equal to p, so it is an identity. 

Now we show that if the state equation p = p(e) holds, the metric 
(12. ip must be static. In fact, from (12. 4p . we know that p = p(r). Using 
the Gj; equation, we obtain 

dtdrOL = 0. 

This implies 

a = f{r) + h{t). 
Replacing t hj t = J e^dt, the metric (12. ip can be rewritten as 

g = -e^/Mrf? + e'^^'-Ur' + r^dO' + sin^ Odi/j''). (2.5) 

Since Gu = Gu (^) = e^°£ if and only if Gu = e'^e, we conclude that 
the fields are static in the sense that they do not depend on i. Thus 
we can assume a = Q;(r) by replacing / by a in (12. 5p . 
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A spherically symmetric star may be assumed to be a ball of radius 
tq centered at r = 0. Then the energy density e and the pressure p 
must satisfy 

e = if r > ro, 
p = if r > tq. 

From now on we assume that the pressure of the star depends on 
r only. Thus the star is static by the above discussion. Let e"^^*^*"^ = 
1 — for some m(r), then 

m(r) = -(r — re~^^). (2-6) 



The G^. equation gives 

da r^p + 2m(r) 



(2.7) 



dr 2r(r — 2m(r)) ' 

and the G* equation leads to 

dm 1 2 
— — = -r e. 
dr 2 

If m(0) 7^ 0, then e^^'-^^ = by (12.61) and the metric g degenerates at 
r = 0. Exclude this case, we may assume 

m(0) = 0. (2.8) 

Finally, in terms of the contracted Einstein equation, we have 

R{g) = e{r) -3p{r). (2.9) 



3. Interior solutions 

In the original investigation on spherically symmetric stars equipped 
with the Tolman metric in [8], the energy density was assumed to be a 
constant. In view of the fact that the energy density should decrease as 
r becomes large, Gu [SBj studied the spherically symmetric dust with 
the energy density given by a (decreasing) step function. In this case 
he found that singularity could only appear away from the origin. The 
complete classification of the spherically symmetric dust was given by 
Hu [36]. In this section, we analyze the results of Section [2] in detail for 
two classes of spherically symmetric stars equipped with the Misner- 
Zapolsky metric [31] and has decreasing energy density functions. 
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3.1. Constant and step function energy densities. We first as- 
sume that 

e{r) = e = constant > for r < ro, (3.1) 
and the mass function is 

m(r) = I K - ""^ f3 2) 

^ ^ I |rQ for r > Tq. 

Then the TOV equation (12. 3p with boundary condition p(ro) = gives 
the following explicit formula for the pressure p{r) inside the star 



p(r) = e , — , ° - for r < rn. (3.3) 

^' 3v/3^^-v/3^^ " 

By Birkhoff theorem [1], the exterior metric must be Schwarzschild 

solution with mass 

D 

S \ / " 



«-t-(^-3rJ*^n^-3fj * (3,4) 
for r > tq. Solving (12. 7p with the continuity condition 

g2a(ro) = 1 _ £!i, (3.5) 

3 

we obtain the interior metric 

2 



5'int 



3 /l_£.2_ 1 /i_£^2] ^^2 



+ (^1 - j ^ dr^ + (d^^ ^ gj^2 ^^^2^ 



(3.6) 



for r < tq. Now it is possible that there exist spherically symmetric 
stars whose energy density e and radius tq satisfy 

3. ^ ^» < 4 p-^> 

Prom (12.91) . we know that the interior metric g^j^^ has singularity at 




n = 3^r2--. (3,8) 

where the pressure function p{r) and the scalar curvature R blow up. 
When To > the singularity appears at 7^ and the origin r = 

is regular. However, for tq = -y/^, we have = and the singularity 
appears at the origin. 
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In [35], Gu used the following step function as energy density to 
study the gravitational collapse: 



Xe 


for 


0<r < A 


fie 


for 


A<r <B 


e 


for 


B < r < tq 





for 


r > ro 



where < A < B < tq, e, X, fi are three positive constants and A > 1, 
< /i < 1. The physical interpretation of the above energy density is 
that the star has a considerable dense shell. 

We shall see that, in certain circumstances, there exists a singularity 
where the energy density is finite but the pressure blows up. 

Set A = ^J^B. Thus we have 

. nA . nB . .B 

^ ' ■ 2 , , / ,.™2 7 _ ^ 2 



Xer dr -\ — / fier dr = — er dr. (3.9) 
2 Jo 2 2 Jo 

Therefore, in the region B < r < r^, the star has mass function 

1 r 

^(^) ~ o e{s)s'^ds 

2 Jo 

.A . i-B . i-r 



Xes ds + - fies ds + - j es ds 

2 ./n 2 ./ /I 2 



B 



B 



A ^ J B 

es^ds 



_ 1 
~ 2 

^ 3 

6 

and the total mass M = |rQ. By Birkhoff theorem P, outside the star, 
the spacetime metric is the Schwarzschild solution (13. 4p . 

Similar to the constant energy density case [8] , plugging the mass for- 
mula in the region B < r < back into the TOV (12. 3p with boundary 
condition p(ro) = 0, we see that the pressure has the same expression 
(13:31). If 



r* = 'i\jrl - — > B, 

there is also a spacetime metric singularity where the energy density is 
finite but the pressure and the scalar curvature blow up. 

3.2. Quadratically decaying energy density functions. In [34j, 

Misner and Zapolsky found an exact solution of the TOV equation 
(12.31) when the energy density function decays quadratically. However, 
their solution does not satisfy the zero pressure condition p(ro) = on 
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the boundary r = tq of the star. Aided by information on the Riccati 
and BernouUi equations, Zhong [37j obtained an exact interior solution 
for the perfect fluid sphere. In both of the above cases, the energy 
density is singular at the origin r = and the physics for such a case is 
somewhat unclear. In this subsection, by modifying the energy density 
function to be finite near the origin, we propose a new model in which 
the energy density is a decreasing function. We shall see that, inside 
the star, there exists a real spacetime singularity and thus the origin is 
not naked. 

Let e and A be two small positive constants such that < £ < |. 
Define the energy density function by 

A for < r < 
= < 4 for ^[^<r < Tq 



A 

for r > ro 
Note that £{r) is a decreasing function and 

lim e{r) = A > — = lim £{r) 



Thus the mass function is 



= < |r for \[^<r < ro 



era 
2 



for r > tq. 

Similar analysis leads to the result that, outside the star (for r > ro), 
the spacetime metric is the Schwarzschild solution 



.e.t--(l-^)*=+(l-?)"<'^= 

+ r^ [de^ + sin^ edip^) 



(3.10) 



with the total mass In the region < r < ro, the metric reads 
^int = -e'^^^'^dt'^ + j^^dr^ + r^dO^ + sin^ Od^^), (3.11) 
where e"'-^-' can be described as follows. Let 

q+ = Vl-e + Vl- 2e, g_ = Vl- s - Vl-2e. (3.12) 

Then 
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Moreover, equation (12 .Sp can be integrated with boundary condition 
p(ro) = to yield the pressure 



e Ir V ^-'^ — Tr ^ 



Pir) = ^ ^ 7=^, (3.13) 



which is valid in the region y ^ < r < tq and is obviously positive. 
The pressure p{r) has finite limit as r — )■ 

The energy density e{r) is discontinuous at r = y^^, but we require 

the pressure be continuous at this point. By contracting the Einstein 
equation (12. 9p . we obtain 

2A 

P_ := lim p{r) = h lim p{r). (3-14) 



In the region < r < i/4f ? the pressure p(r) satisfies the equation 



dp dr 



(A + p(r))(A + 3p(r)) 2(Ar2 - 3) 
with boundary condition 



(3.15) 



I 2A ^ 



;f ) V - - 



l-2e -2, 

rn 



^ 9 ' „ r „ n „ -2 / 1-2^-1 



3£ 

A 



As in the case of constant energy density, one can obtain an explicit 
solution of equation (I3.15P with singularity at 



9(3-3£)(P_ + A)2 



A(3P_ + A)^ 



(3.16) 



Curvature tensors are finite but discontinuous at r = y y . So the 

point r = is not a spacetime singularity. 

Note that the density function in reference [37] is of the form p = 
^Q^Q^2 ( PTl equation (4)]), where c and G denote the speed of light and 
Newton's constant respectively. Thus, the coefficient of ^ is fixed and 
the only free parameter is the star radius tq. In our context, both e 
and ro are free parameters. If we set the parameter e = |, then the 
exterior Schwarzschild solution (I3.10p with the total mass ^ = ^ can 
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be glued to the interior solution f l3.1ip at the boundary r = ro of the 
star, i.e. 

g = + ^ + {de^ + sin^ ed^^) . 

This recovers (27) and (28), the particular case considered in [37] . 
4. Quantising interior solutions 

In this section we quantise the metrics (13. 6p (see Section [3]) following 
the approach to noncommutative general relativity developed in [TTl 
[T2| [T3] . The purpose of our investigation is to examine possible effects 
of space-time noncommutativity by considering simple models. 

Note that the interior metrics (13.61) and (13.111) for spherically sym- 
metric stars can both be written as 

gmi = -a^{r)dt^ + l?{r)dr'^ + r'^{de'^ + siYi^ed-ip'^), r < ro, (4.1) 

for appropriate functions a(r) and h{r) where 

a(r,) = 0, a'ir,) ^ 0, a"(r,) ^ 0, b{r,) ^ 0, b'{r,) ^ (4.2) 

for some r^. < Tq. 

To carry out the deformation quantisation (in the sense of [HI [T21 
[13] ) of the the metric (14.11) . we need first to specify a Moyal algebra. 
Denote = t, = r, x'^ = 6 and x^ = (p. We deform the algebra of 
functions in these variables by imposing the Moyal product 

(/ * g){x) = f{x) expt^"'^^^ g{x) 
with the following anti-symmetric matrix 



/ 








o\ 























1 







-1 


0/ 



(4.3) 



where h is the deformation parameter, which may be regarded as re- 
lated to the Planck constant. 

Some comments are in order. Noncommutativity of time coordinate 
results in violation of unitarity and causality for quantum field theories 
defined on flat spacetimes as well as spacetimes with compact spatial 
submanifolds [30] (also see [3l]; for a discussion in the context of string 
theory, see [32]). In order to retain basic principles of quantum physics 
such as unitarity and causality, we have to keep the time coordinate 
commutative, and the Ansatz (14.31) enforces this. Note also that the 
Ansatz (14.31) leads to the simplest possible model of space-time non- 
commutativity. 
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Following [llj, we choose the following embedding 

X^ = a{r)smt, X^ = a{r)cost, = f{r), 

4 5 6 (4-4) 

X = r sin ^ cos X = r sin ^ sin X =rcos6. 



where /(r) is related to the functions a(r) and b{r) through the differ- 
ential equation 

{fr+i = {ar+b'. 

Here the rationale is much the same as in the classical theory of 
surfaces, which can be treated as two-dimensional geometries embedded 
in the three-dimensional Euclidean space. The surfaces themselves are 
determined by the embeddings, and their differential geometry can be 
studied by analyzing the embeddings in a very concrete manner using 
elementary techniques (see, e.g., |33] for an introduction to the theory 
of surfaces) . 

The embedding f l4.4p determines the noncommutative geometry of 
the space-time, which we now analyze. The quantum deformation of 
the metric (14. ip is defined by 

g^, = - d^X' * d,X' - d^X^ * d^X^ + d^X' * d,X^ 

+ d^X* * d,X^ + d^X'' * d,X'> + 9^X6 * '^^ 

Note that in the classical commutative limit with /i = 0, (14. 5 p reduces 
to g^,u = -Ei=idf.X'd,X' + E,=3'9/.^^'5.^^ which indeed recovers 
the metric (14.10 as one can verify. 

Lengthy computations yield the following result for the noncommu- 
tative metric (14.50 : 

goo = - a'^ir), 

gOl =glO = g02 = g20 = g03 = gSO = 0, 

gii =b'^{r) + (sin^ 6 — cos^ 6) sinh^ 
gi2 =g2i = 2r sin 9 cos 6 sinh^ h, 
gi3 = ^ gai = ~2r sin 6 cos 6 sinh h cosh h, 
g22 =r^ [l - (sin^ 6 - cos^ 6) sinh^ h] , 
g23 = — g32 = r'^ (sin^ 9 — cos^ 9) sinh h cosh h, 
g33 =r^ [sin^ 9 + (sin^ 9 — cos^ 9^ sinh^ . 

Given the general theory of |JJJ, the computation of the connection 
and Riemaniann curvature of the noncommutative metric (14.60 is in 
principle straightforward though very laborious. We shall not spell out 



(4.6) 
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the complete results here. However, we shall need the noncommutative 
scalar curvature, which can be expressed in the form 

2^2^3 ra{r)A^ a{r)A ' ^ ' ' 

where A, Aq, Ai and A2 are given in the Appendix. 

For the sake of being concrete, we shall only discuss the noncom- 
mutative singularity for quantum analogue of the metric (13. 6p . As we 
have seen in Section |31 singularity occurs essentially in the region with 
constant energy density. Similarly in the present setting, if the en- 
ergy density e and radius tq of a spherically symmetric satisfy (I3.7jl . 
noncommutative singularity also occurs for the metric (14. ip at 



Indeed, in this case. 



a(r) = - 4 / 1 — -rl — - 1 / 1 — -r^, 
^ ^ 2V 3 ° 2V 3 



h{r) = {l-^r 

we obtain 



_ 1 

e 2\ 2 



(4.8) 



R = 5i - §, (4.9) 

where 5, Bi and B2 are given in the Appendix. 

Note that -B|r=r* = 0. Inspecting the formulae for Bi\y.=^^ and 
i?2|r=r, in the Appendix, we see that Bi\^=y.^ is finite and 52|r=r, 7^ 0. 
Thus we have 

R(n) = 00, (4.10) 

and hence r = r^, is a singularity in the noncommutative spacetime. 

Remark 4.1. If r^e > |, then r = is a regular point of the quantum 
deformed metric for 1^3. 6\) . In this case, the noncommutative scalar 
curvature 

e fcosh 2h - 2^9- -irle + 3) sech% 

R= ^ , ^ 

9 - 'irle - 1 

is finite. When tq = \/^, the origin r = is singular. 
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5. NONCOMMUTATIVE GRAVITATIONAL COLLAPSE 

In this section, we quantise the dust solutions [8] and study noncom- 
mutative gravitational collapse. While our method for quantisation is 
much the same as in Section HI a new feature is that the quantised dust 
solutions have an explicit time dependence and their time evolutions 
are thus clear. 

By replacing r by t, and R by r, the Tolman metric studied in [8] 
can be written as 

ds^ = -de + (1 - ctf/'' [dr^ + r\de^ + sin^ ed<p'')\ (5.1) 

1 _3 

with c = 3rg i?^ ^ , where tq is the gravitational radius and Rh is the ra- 
dius of the star or some micro object. This spacetime can be embedded 
into a 5-dimensional flat Minkowski spacetime via 

%l-ctf'^ (r^ 
X^=^^^+(^ + l)(l-ct)n 

. 9(1 -ct)^/3 .^2 

=(1 - ct)2/Vcos0sin^, 

=(1 - ct)2/3r sin ^ sin 0, 

=(1 -ct)2/3r cos ^. 

As in Section HJ we deform the algebra of functions in the variables 
r, t, TT and ^ into a Moyal algebra ^ defined by the anti-symmetric ma- 
trix (14. Sp . Now we consider the noncommutative geometry embedded 
in Jc' by (15. 2p . The noncommutative metric of the embedded noncom- 
mutative geometry (defined in the standard way [11]) yields a quantum 
deformation of the metric (15. ip : 

'3 



g/.. = - d^X^ * d,X' + d^X^ * d,X^ + d^X^ * d,X' 

+ d^X^ * d,X^ + d^X'' * d,X\ 

which can be computed explicitly. We have 

4c^r^ cos 29 sinh^ h 
" 9(1 - ct)2/3 ^' 

2 

gi2 =g2i = -cr(l — ct)^^'^ cos 20 sinh^ 
3 

4 

gi3 =g3i = — -cr^(l — ct)^" cos 6'sin 6'sinh^ /i, 

gi4 = - g4i = icr^(l - ct)^/^sin20sinh2/;., 
3 



(5.3) 
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g22 =(1 - ct)^/^ (l - COS 26 sinh^ h) , 
g23 =g32 = ^(1 - ct)'^/^sin26'sinh^ h, 
g24 = — g42 = — 2r(l — cty^^ cos 6 cosh h sin ^ sinh h, 



g33 =r2(l - ct)^/^ (cos 2^ sinh^ + l) , 

g34 = - g43 = -^r^l - ct f^ cos29smh2h, 

g44 = - -r^(l - ct)^/^(cos20cosh2/i - 1). 



The noncommutative scalar curvature is given by 

_ 4c2 cosh^ /z Ci 
- (1 _ct)4/3C3 



(5.4) 



where C and Ci are given in the Appendix. 

We shall consider aspects of the noncommutative spacetime by ex- 
amining the behaviour of the scalar curvature as time increases by 
following an approach adopted in [5] in the classical context. When 
time approaches values where R — )■ oo, the radius of the stellar ob- 
ject reduces to zero, and this is an indication of gravitational collapse 
[5]. Obviously this only provides a snapshot, nevertheless, it enables 
us to gain some understanding of gravitational collapse in the non- 
commutative setting. A full treatment of the time evolution of stellar 
objects ending at gravitational collapses in noncommutative geometry 
will be given in a future publication. 

Let us regard h as a real number and make the (physically realistic) 
assumption that h is positive but close to zero. Now if t is significantly 
smaller than ^ compared to h, that is, ^ — t ^ h, both the noncom- 
mutative metric and noncommutative scalar curvature R are finite, 
and there is non-singularity in the noncommutative spacetime. Thus 
the stellar object described by the noncommutative geometry behaves 
much the same as the corresponding classical object. 

When t = t^. := -, we have R|t =i = cxd and the radius of the stellar 

^ * c 

object reduces to zero. This is the time when gravitational collapse 
happens in the usual classical setting. 

However, in the noncommutative case, singularities of the scalar cur- 
vature already appear before t*. Indeed, when time reaches 




c 



1 



1 




(2 cos 2^ + cosh 2/1 + 3)^/2 



cosh h 
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C vanishes and Ci/(1 — ct)^/^ is finite of order in h. Tlius tlie scalar 
curvature tends to infinity for all t(r, 9) and tlie noncommutative space- 
time becomes singular. 

Tliis indicates tfiat in tfie noncommutative setting, gravitational col- 
lapse happens within a certain range of time because of the quantum 
effects captured by the non-commutativity of spacetime. This is fully 
consistent with usual expectations of quantum mechanics. However, 
effect of non-commutativity only starts to appear at third order of h. 



Appendix A. Some formulae used in the main text 

This appendix spells out some lengthy formulae which have been 
used in the main body of the paper. The quantities A, Aq, Ai and A2 
used in f l4.7p are given by 



Z cos 2^ + cosh2/i + 3) sinh^/i. 



A =2(^cos2^sinh2/i - ljb{rf 

Ao = - 8 (^3 cos 26 sinh^ 2h + 2 cosh 2h + 2 cosh 4/i j h{rf 

+ 2 ( - 8 cos 4^ sinh^ /i + 92 cosh 2h - 19(cosh 4/i + 3)) cosh^ h 
+ 4 cos 2^(7 - 4 cosh 2h) sinh^ 2h\ hirf 



+ 2r 



(8 cos AO sinh^ h — 28 cosh 2h + cosh Ah — 5) cosh^ h 



+ 2 cos 26(3 cosh 2h - 1) sinh^ 2h b'{r)b{r 



-A 



4( cos 2^(3 cosh 2/i - 23) - 2 cos 4^) cosh^ /i sinh^ /i 



(-14cosh2/i + cosh4/i + 9) sinh^2/i 6(r 



cos 2^(3 cosh Ah + 13) sinh^ 2h 



- 16(cos A9 - cosh 2h + 2) cosh^ h sinh^ h h'{r)h{r) 
+ 16-7 cos 29 - cos A9 + (3 cos 29 + 4) cosh 2h - h 

Ai =8 cos 29 cosh^ h sinh^ h + (3 cosh 2h + 1 j sinh^ 2h 
+ 2 (6 cosh 2h + 2) cosh^ h + cos 29 sinh^ 2h h{rf 

— 8r^cos2^sinh^ /i — ij b{r)b'{r), 
Ao =a( cos2^sinh2/i - iV 



cosh^ h sinh^ h. 
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The following quantities appeared in equation (14. 9p : 



B = ~3 + er^ + 3^(3 - r%)(3 - r^e) 

e 



sinh h — 6 



2r% cos 26 + (r% - 3) (cosh 2h + 3) 

- IGr^e^ cos 49 sinh^ h - 2r^e {23er^ + 9) 
+ {er^ (68^2 + 2l) + 405) cosh2/i + 2 {er^ (9 - Ur'^e) + 8l) cosh4/i 
+ (r^e - 3) (4r^£ - 9) cosh Qh + 270 



cosh h 



+ 2r^e cos 2^(3 cosh 2h - 7) ((r^e - 3) cosh 2h - 3 - er'^) sinh^ 2/i| , 



5, = - 



1 



2r% cos 29 + (r% - 3) (cosh 2h + 3) 



sinh h — Q 



e (r^e — 3) cosh^ h 



£ cos 29 sinh h — r e 



+ 4 (er^ + 9) cosh 2/i + (9 - 3r2e) cosh 4/i + 27 | . 



Evaluating Bi and i?2 at r 



3^/^o — ^, we obtain 



5i| 



6 {3rl£ - 8) cos 20 + 9 (r^e - 3) (cosh 2h + 3) 



X <^ 9 cosh h 



sinh /i — 6 



[<ol2e\l - 32A3erl + 4341) cosh2/i 

- 2 (3e (39r2£ - 21l) + 847) cosh Ah 
+ 3 (r^e - 3) {I2rle - 35) cosh Qh 

- 2 (^8 (8 - 3r2e)^ cos40 sinh^ /i + 3rle (69r^e - 365) + 1433^ ^ 



+ 18 sinh^ 2h 



{3rlE - 8) cos 29 (3 cosh 2/i - 7) 



{-^erl + 3 [rle - 3) cosh 2/i + 7) } 
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2 ^ Q\ r^t-.c^'U, 2 



3e {tqE — 3) cosh h 

-D2 r=r. 



( 2 (Sr^e - 8) cos 29 + 3 {r^e - 3) (cosh 2h + 3) sinh^ h - 2^ 
X 8 (8 - 3r^£) cos 26 sinh^ /i - 3rge 

+ 4 (3r^e - 5) cosh 2h-9 {rle - 3) cosh Ah + 17 . 

These formulae are used in the derivation of fl4.10p . 

The quantities C and Ci in (15 ■4p in the main body of the paper are 
given by the following formulae: 

C =9(1 - ct)^/^ cosh^ h - 2cV2(2 cos 29 + cosh 2h + 3) sinh^ h, 

Ci = - 243(1 - ct)^/^ cosh*^ h + 486(1 - ct)^/^ cosh^ h 

- 18cV2(l - ctf/^{2 cos 2^-3 cosh 2h - 1) sinh^ h cosh^ h 

- QcV^ (1 - cty/^ (52 cosh 2/1 + 3 cosh Ah 

+ cos 2^(28 cosh 2h + cosh 4/i - 13) + 9^ sinh^ h cosh^ h 
+ Ac^r^ sinh^ h (a cos 2^(cosh 2h + 15) sinh^ /i 
+ 2 cos 4^(cosh 2/i - 3) + 38 cosh 2h + 3 cosh 4/i - 13^ 



2 
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